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Abstract 

In this paper, we calculate the branching ratios and CP-violating asymmetries for B s — > 
p^K^, B s — > p°K and B s — > ujK decays in the perturbative QCD (pQCD) factorization 
approach. The theoretical predictions for the CP averaged branching ratios of the considered 
decays are: Br(B s -> p ± K T ) ps 24.7 x 10" 6 , Br(B s p°K°) ps 1.2 x 10~ 7 and Br(B s 
wF) « 1.7 x 10~ 7 ; and we also predict very large direct CP-violating asymmetries for the 
latter two decay modes: A d ^ P {p ± K^) » -12%, A% P (p°K°) « -92%, A^ P {co K°) ps 81%, 
A^|?(p K ) ps -36%, and ps -40%, which can be tested in the forthcoming LHC- 

b experiments. 

PACS numbers: 13.25.Hw, 12.38.Bx, 14.40.Nd 
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I. INTRODUCTION 



The charmless B meson decay is a good place to test the Standard Model (SM), study 
CP violation and look for signal or evidence of possible new physics beyond the SM. 
Since 1999, many such decay modes have been observed in the B factory experiments. In 
the forthcoming Large Hadron Collider beauty (LHC-b) experiments, a large number of 
heavier B s and B c mesons together with light B Ut d mesons will be produced The study 
about the charmless decays of B s meson therefore is therefore becoming more interesting 
then ever before 0. 

By employing thegeneralized facterization approach 0, 0, 0] or the QCD factorization 
(QCDF) approach p, about 40 B s — > h\h 2 (hi stand for light pseudo-scalar or vector 
mesons ) decay modes have been studied in the framework of SM Fl @Jp, Ell or m some 
new physics models beyond the SM In Refs. 0, 0, Q, fl5l lid], the branching 



ratios and CP violating asymmetries of B s — > n + n~,7rp, ttK, p(uj)K* and nr]^'' decays 
have been calculated by employing the perturbative QCD (pQCD) factorization approach 
fl7L [lH . Motivated by the expected successes in LHC-b experiments and other hadronic 
B meson experiments, we here continue to investigate more charmless B s decays in pQCD 
factorization approach. 

In this paper, we will study the B s — > p ± K T , p°K and uK decays in the pQCD 
approach. In principle, the physics for the B s two-body hadronic decays is very similar 
to that for the B d meson except that the spectator d quark is replaced by the s quark. 
Theoretically, the three decays have been studied before in the naive or generalized fac- 
torization approach, as well as in the QCD factorization approach |gL [To| : 

For B s — > p{u))K decays, the B s meson is heavy, setting at rest and decaying into 
two light mesons (i.e. p(u) and K ) with large momenta. Therefore the light final state 
mesons are moving very fast in the rest frame of B s meson. In this case, the short 
distance hard process dominates the decay amplitude. We assume that the soft final 
state interaction is not important for such decays, since there is not enough time for light 
mesons to exchange soft gluons. Therefore, it makes the pQCD reliable in calculating 
the B s — > p{uj)K decays. With the Sudakov resummation, we can include the leading 
double logarithms for all loop diagrams, in association with the soft contribution. Unlike 
the usual factorization approach, the hard part of the pQCD approach consists of six 
quarks rather than four. We thus call it six-quark operators or six-quark effective theory. 
Applying the six-quark effective theory to B s meson decays, we need meson wave functions 
for the hadronization of quarks into mesons. All the collinear dynamics are included in 
the meson wave functions. 

This paper is organized as follows. In Sec. |H1 we give a brief review for the pQCD fac- 
torization approach. In Sec. IHIl we calculate analytically the related Feynman diagrams 
and present the decay amplitudes for the studied decay modes. In Sec. IIV1 we show the 
numerical results for the branching ratios and CP asymmetries of B s — > p(ui)K decays 
and comparing them with the results obtained in the QCDF approach. The summary 
and some discussions are included in the final section. 
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II. THEORETICAL FRAMEWORK 



The three scale pQCD factorization approach has been developed and applied in the 
non-leptonic P( s ) meson decays for some time. In this approach, the decay amplitude 
is separated into soft ($), hard(H), and harder (C) dynamics characterized by different 
energy scales (t, m&, Mw) . It is conceptually written as the convolution: 

A(B (S) -> M X M 2 ) ~ Jd%d%d% Tr [c{t)$ B{s) {k 1 )$M 1 {k2)®MAh)H{k x MM>t) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the 
trace over Dirac and color indices. The harder dynamic involves the four quark operators 
described by the Wilson coefficient C(t). It results from the radiative corrections to 
the four quark operators at short distance. In the above convolution, C(t) includes the 
harder dynamics at larger scale than Mg (s) scale and describes the evolution of local 4- 



Fermi operators from myy (the W boson mass) down to t ~ 0{^J KMb {s) ) scale, where 

A = Mb, — rrii,. The function H(k\, k 2 , ks,t) describes the four quark operator and the 
spectator quark connected by a hard gluon whose scale is at the order of Mg (s) ,so this hard 
part H can be perturbatively calculated. The hard and harder dynamics together make 
an effective six-quark interaction. The soft dynamic is factorized into the meson wave 
function &m ,which describes hadronization of the quark and anti-quark pair into the 
meson M. While the function H depends on the processes considered, the wave function 
$a/ is independent of the specific processes. Using the wave functions determined from 
other well measured processes, one can make quantitative predictions here. 

For the B s meson decays, since the b quark is rather heavy we consider the B s meson 
at rest for simplicity. It is convenient to use light-cone coordinate (p + ,p~ , p T ) to describe 
the meson's momenta. The B s meson and the two final state meson momenta can be 
written as 

^ S = ^(1,1,0 T ), P p(w) = ^(l,rJ (w) ,O r ), P fc = ^(0,l-rJ H ,O T ), (2) 

respectively, where r p = m p /M Bs ; the light pseudoscalar meson masses have been ne- 
glected. 

For the B s — > p{oj)K decays considered here, only the p(u) meson's longitudinal part 
contributes to the decays, its polar vector is = J^®^ - (1, —r 2 ^, T ). Putting the light 

(anti-) quark momenta in B s , p(u) and k mesons as fci, and ks, respectively, we can 
choose 

h = (xiP+,0,k 1T ), k 2 = (x 2 P+,0,k 2T ), k 3 = (0,x 3 P 3 -,k3 T ). (3) 

Unlike the QCD factorization approach,we don't neglect the transverse momentum kx in 
the above expressions, by which to avoid the endpoint singularity. Then, the integration 
over k{ , k~2~ , and k 3 in eq.((H) will lead to 



A(B S — > p{oj)k) ~ J dxidx 2 dx 3 bidbib 2 db2b 3 dbs 

Tr ^(^^(x!,^)^)^,^)^^,^)^^,^^)^^^-^] , (4) 
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where b{ is the conjugate space coordinate of kiT, and t is the largest energy scale in 
function H(xi,bi,t), as a function in terms of Xi and cV The large logarithms (hxmw/t) 
coming from QCD radiative corrections to four quark operators are included in the Wilson 
coefficients C(t). The large double logarithms (ln 2 x«) on the longitudinal direction are 
summed by the threshold resummation and they lead to S t (xi) which smears the 

end-point singularities on x«. The last term, e~ s ^\ is the Sudakov form factor resulting 
from overlap of soft and collinear divergences, which suppresses the soft dynamics effec- 
tively [i^. Thus it makes the perturbative calculation of the hard part H applicable at 
intermediate scale, i.e., Mg s scale. We will calculate analytically the function H(xi,bi,t) 
for B s — > p{oj)K decays in the first order in a s expansion and give the convoluted ampli- 
tudes in next section. 



Wilson Coefficients 



It is well known that the low-energy effective Hamiltonian is the basic tool to calculate 
the branching ratios and Acp of B meson decays. For B s — > p{uj)K decays, the related 
weak effective Hamiltonian H e ff can be written as jijj 



n 



eff 



G p 
71 
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v ub v: d (c a (/i)o?o*) + c 2 (p)o u M) - v tb v; d £ c^) o^) 



i=3 



(5) 



We specify below the operators in H. e ff for b — > d transition: 



0\ 

7 
9 



\dyfLbc 
Ida^Lbc 



Eg' e q ,q'^ IM Rq'p , 
Eg' e q/ tf /3 ^Lq J p , 



01 

O e 
8 
O w 



dyfLbp ■ Eg' q'/3lv L( la > 
dcffLbp ■ Eg' tfplfiRq'a , 
\d a ^Lbfi ■ E 9 / e ql q'^^Rq' a 
\d a ^Lbp ■ Eg' z q >q'f3lnLq' a 



(6) 



where a and /3 are the SU(3) color indices; L and i? are the left- and right-handed 
projection operators with L = (1 — 75), R = (1 + 75). The sum over q' runs over the 
quark fields that are active at the scale \i = 0(m b ), i.e., (q'e{u, d, s, c, b}). The pQCD 
approach works well for the leading twist approximation and leading double logarithm 
summation. For the Wilson coefficients Cj(/i) (i = 1, . . . , 10), we will also use the leading 
order (LO) expressions, although the next-to-leading order calculations already exist in 
the literature |2j]|. This is the consistent way to cancel the explicit /x dependence in the 
theoretical formulae. 

For the renormalization group evolution of the Wilson coefficients from higher scale to 
lower scale, we use the formulae as given in Ref. [2^ directly. At the high mw scale, the 
leading order Wilson coefficients Ci(Mw) are simple and can be found easily in Ref. [2JJ . 
In pQCD approach, the scale t may be larger or smaller than the m b scale. For the case 
of m b < t < mw, we evaluate the Wilson coefficients at t scale using leading logarithm 



running equations, as given in Eq.(Cl) of Ref. 2^1 . For the case of t < m b , we then 
evaluate the Wilson coefficients at t scale by using Ci{m b ) as input and the formulae 
given in Appendix D of Ref. 2^ . 
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B. Wave Functions 



In the resummation procedures, the B s meson is treated as a heavy-light system. In 
general, the .B^rneson light-cone matrix element can be decomposed as [a. l23l| 



' d'z ikl . z 



, e^(0\b a (0)s^z)\B(p Bs )) 



s s (ki) - — ^0 B .(ki) 



(7) 

/3a 



where n = (1,0, T ), and v = (0, 1, T ) are the unit vectors pointing to the plus and 
minus directions, respectively. From the above equation, one can see that there are two 
Lorentz structures in the B s meson distribution amplitudes. They obey to the following 
normalization conditions 

In general, one should consider these two Lorentz structures in calculations of B s meson 
decays. However, it can be argued that the contribution of <pB a is numerically small, thus 
its contribution can be numerically neglected. Using this approximation, we can reduce 
one input parameter in our calculation. Therefore, we only consider the contribution of 
Lorentz structure 

$Bs = -^W^ Bs + M ^)^^(ki)- (9) 

In the next section, we will see that the hard part is always independent of one of the 
kf and/or k± , if we make approximations shown in next section. The B s meson wave 
function is then the function of variable k^ (or kf) and kj~, 

<j> Bs (K,kt) = J dkMk+,K,kt). (io) 

The wave function for the pseudoscalar meson K are given as: 

$ K (P, x, C) = -^=75 {Mi(x) + m^ p K (x) + Cm^{fi - v ■ n)f K (x) } (11) 
V Zl\ c 

where P and x are the momentum and the momentum fraction of K meson, respectively. 
We assumed here that the wave function of K meson is the same as the wave function 
of 7r meson. The parameter ( is either +1 or —1 depending on the assignment of the 
momentum fraction x. 

For B — ► pK decay, the p meson is lon gitu dinally polarized. We only consider its wave 
function in longitudinal polarization jiil 125^ . 

1 rl 



< p-(Pe L )\d a (z)u p (0)\0 >= -= / dxe™ p - z y [j p( f (x) + m p <j> p {x)] + m p <f> s p (x)}(12) 

V ZI ^c Jo 
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The second term in above equation is the leading twist wave function (twist-2), while the 
first and third terms are sub-leading twist (twist-3) wave functions. For B — > ujK decay, 
we have similar expression as Eq. (|12|). 

The transverse momentum k x is usually conveniently converted to the 6 parameter by 
Fourier transformation. The initial conditions of function <pi(x), i = B s ,p,u,k, are of 
non-perturbative origin, satisfying the normalization relation 

f 1 1 

J <f>i{x, b = 0)dx = -^j^fi , ( 13 ) 
with fi the meson decay constants. 



III. PERTURB ATIVE CALCULATIONS 

In the previous section we have discussed the wave functions and Wilson coefficients of 
the amplitude in eq.(fTj). In this section, we will calculate the hard part H(t). This part 
involves the four quark operators and the necessary hard gluon connecting the four quark 
operator and the spectator quark. We will show the whole amplitude for each diagram 
including wave functions. Similar to the B s — > nK the total 8 lowest order diagrams 



contributing to the B s — ► pK decays, are illustrated in Figure 1. We first calculate the 
usual factorizable diagrams (a) and (b). Operators Oi, O2, 3 , O4, O9, and 0\q are 
(V — A)(V — A) currents, the sum of their amplitudes is given as 



FeK = -8nC F f p M^ g / dx ± dx 3 \ &id&iM&3 0B.(zi, &i) 

Jo Jo 

• { [(2 - x 3 )(j>£(x 3 , 63) + (2x 3 - l)r K (^(x 3 , b 3 ) - f K {x 3 , 63))] 
a s (tl)h e (xi, 1 - x 3 , h, 63) exp[-S a b(tl)] 

+2r k (j)^(x 3 ,b 3 )a s (tl)h e (l -x 3 ,x u 63,61) exp[-^ a6 (^)]} , (14) 

where Ci? = 4/3isa color factor. The function h l e , the energy scales tf e and the Sudakov 
factors S a b are displayed in the Appendix. In the above equation, we do not include the 
Wilson coefficients of the corresponding operators, which are process dependent. They 
will be shown later in this section for different decay channels. The diagrams Fig. 1(a) 
and 1(b) are also the diagrams for the B s — > K form factor F Q S ~* K . Therefore we can 
extract F Bs ^ K from Eq. (fHjl . 

The operators O5, Oq, O7, and 0$ have a structure of (V — A) (V + A). In some decay 
channels, some of these operators contribute to the decay amplitude in a factorizable way. 
Since only the axial- vector part of (V + A) current contribute to the pseudo-sealer meson 
production, 

(p\V - A\B)(K\V + A\Q) = -(p\V - A\B)(K\V - A\0), (15) 

the result of these (V — A)(V + A) operators is opposite to Eq. (JHJ), i.e., 

F P eK = -F eK . (16) 

In some other cases, we need to do Fierz transformation for these operators to get right 
color structure for factorization to work. In this case, we get (S + P)(S — P) operators 
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B x 




(a) 



K B s 




(b) 



A' 



B x 



(c) 



K B s 



(d) 



K 



B x 




B x 




B x 




FIG. 1: Diagrams contributing to the B s — » pA and wif decays (diagram (a) and (b) contribute 



to the — > K form factor F, 



o 



from (V — A)(V + A) ones. For these (S + P)(S — P) operators, Fig. 1(a) and 1(b) will 
give 



7 p 2 
eK 







(17) 



For the non-factorizable diagram 1(c) and 1(d), all three meson wave functions are 
involved. The integration of b 3 can be performed using 5 function 5(b 3 — 61), leaving only 
integration of b\ and b 2 . For the (V — A)(V — A) operators, the result is: 

32 f 1 f 00 

M eK = —irC F Mg 3 / dx\dx 2 dx 3 / b 1 db 1 b 2 db 2 (J)B s (xi,bi)(j)p(x 2 ,b 2 ) 
V6 3 Jo Jo 

• {(1 - x 3 ) [(t)£(x 3 , h) + 2r k( f)l(x 3 , h)] 

■ oc s (t f )hf(xi,x 2 , 1 - x 3 , h, b 2 ) exp[-S cd (t f )} } . (18) 
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For the (V — A)(V + A) operators, the result is: 

64 f 1 f°° 

M^k = —=ttC f G f M b r p \ dx x dx 2 dx 3 \ b 1 db 1 b 2 db 2 (f)B s ( x i^i) 
v6 Jo Jo 

■ { [x 2 <Pk{x^ 61) (<t> s p {x 2 , b 2 ) - </) P {x 2 , b 2 )) + r K ((1 + x 2 - rr 3 ) 

(0x(^3, &l)0p(»2, 62) - 0x(^3, &l)$(&2, 62)) + (1 - X 2 - X 3 ) 

(0x(^3, &i)0*(z 2 , 62) - 0x0^3, &i)^(&2, 62)))] 

a s (t f )h f (x u x 2 , 1 - x 3 , 61, 6 2 ) exp[-S' C(i (t / )] } . (19) 

For the (S* — P)(S + -P) operators, the result is: 

M% = -M eK . (20) 

For the non-factorizable annihilation diagram 1(e) and 1(f), again all three wave func- 
tions are involved. Here we have two kinds of contributions. M a k is the contribution 
containing operator type (V — A)(V — A), while is the contribution containing op- 
erator type (V -A)(V + A): 

32 f 1 f°° 

M aK = —nC F M B I dx 1 dx 2 dx 3 / 6id6i6 2 d6 2 <Pb s (xi, h) 
v6 Jo Jo 

■ { [x 2 (j) p (x 2 , b 2 )(pj ( (x 3 , 62) + r p r K ((x 2 - x 3 ) 

■ {<Pk(^3, b 2 )4> P {x 2 , b 2 ) + f K {x 3 , b 2 )<j) s p {x 2 , 6 2 )) + {2 + x2 + x3) 

■ <Pk(%3, b 2 )(f) s p (x 2 , b 2 ) + (-2 + x 2 + x 3 )(f)l(x 3 , b 2 )^ p (x 2 , b 2 ))] 

■ a a (tf)h 2 f (x 1 ,x 2 , x 3 , 61, 6 2 ) exp[-S e f(tf)] 

+ [~x 3 (pp(x 2 , 62)0^(^3, 62) + r p r K ((x 2 - x 3 ) (0£(z 3 , b 2 )(j) p (x 2 , b 2 ) + (j) T K {x 3) b 2 ) 

■ <Pp{x 2 , b 2 )) - (x 2 + x 3 ) (0£(x 3 , b 2 )<f) s p {x 2 , b 2 ) + 4>k{x 3 , 6 2 )0* (x 2 , b 2 )))] 

■ a s (t 1 f )h 1 f (x 1 ,x 2 , x 3 , 61,62) exp[S e f(tf)]} , (21) 

32 f 1 f°° 

M ak = —j=^C F M% I dx x dx 2 dx 3 \ bxdb x b 2 db 2 (f>B e (xi, h) 
V6 3 Jo Jo 

■ { [r p {x 2 - 2)<f ) £{x 3 , b 2 ) (0* (x 2 , b 2 ) + <j) s p {x 2 , b 2 )) 
+r K (2 - x 3 )(p p (x 2 , b 2 ) (0£(a?3, h) + <Pk(x 3 , 62))] 



a s (t 2 f )h 2 f (xi, x 2 , x 3 , 61, b 2 ) exp[-S ef (tf)] 



+ [r p (-x 2 )(f)£(x 3 ,b 2 ) (0*(rr 2 ,6 2 ) + 0*(x 2 ,6 2 )) 
+r K x 3 <pp(x 2 , b 2 ) (<Pk( x 3, h) + 4> T K {x 3 , b 2 ))] 

■ a s (t 1 f )h 1 f (x 1 ,x 2 , x 3 , 61,62) exp[-S e f(if)]} , (22) 

where tk = /Mb s and r p = m p /MB s - 

The factorizable annihilation diagram 1(g) and 1(h) involve only p and K wave func- 
tions. There are also two kinds of decay amplitudes for these two diagrams. F aK is for 
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(V — A)(V — A) type operators, and is for (S — P)(S + P) type operators, 

P 1 POO 

FaK = 8irC F f Bs M Bs / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

■ { [x 3 (j) p (x 2 , b 2 )(f)^(x 3 , b 3 ) + 2r p r k (f) s p (x 2 , b 2 ) ((1 + x 3 )(j)^(x 3 , b 3 ) 
(x 3 - \)4> T K (x 3 , b 3 ))] a s (t 3 e )h a (x 2 , x 3 , b 2 , b 3 ) exp[-S gh (t 3 e )} 

- [x 2 (f)p{x 2 , b 2 )(pK(x 3 , b 3 ) + 2r p r k (f)^(x 3 , b 3 ) ((1 + x 2 )(f> s p (x 2 , b 2 ) 

(x 2 - 1)0* (^2, b 2 ))] a s (te)h a (x 3 ,x 2 , 6 3 ,6 2 ) exj)[-S gh {ti)]} , (23) 



f p 2 



— 167tCf/s 3 M^ / dx 2 dx 



b 2 db 2 b 3 db 3 

■ { [2r p (p p (x 2 , b 2 )(j) K (x 3 , 63) + r k x 3 (p p (x 2 , b 2 ) (<fr K (x 3 , 6 3 ) 
-<j) T K {x 3 , 63))] a s (t 3 e )h a (x 2 , x 3 , 6 2 , 63) exp[-Sg h (t 3 e )] 

+ [x2rp(p K (x 3 , 63) (4> s p (x 2 , 6 2 ) - 0* (x 2 , 6 2 )) + 2r k (p p (x 2 , b 2 ) 

■ (p K (x 3 ,b 3 )a s (t 4 e )h a (x 3 , x 2 ,b 3 ,b 2 ) exp[-S gh (tj)]} . 



(24) 



In the above equations, we have assumed that X\ « x 2 ,x 3 . Since the light quark 
momentum fraction x% in B meson is peaked at the small X\ region, this is not a bad 
approximation. The numerical results also show that this approximation makes very little 
difference in the final result. After using this approximation, all the diagrams are functions 
of hi = x x M B J^[2 of B s meson only, independent of the variable of kf. Therefore the 
integration of eq. ffTU|) is performed safely. 

Combining the contributions from different diagrams, the total decay amplitude for 
B s — > p + K~ decay can be written as: 



M{p + K~ 



where £ u = V* b V ud , & 



-C 3 + C4 + -Cg + C10 

6 (C 3 + C 9 )] + MS Kt (C 5 + C r )] 



+M aK 
+F a x 
+F a P K 

v; b v td 



-6 
-6 



C3 — — Cq 

6 2 

\c 3 + C 4 • 



+ M a P ^ 



C 5 - -C 7 



*C 5 
3 



— Cq — — Cin 

6 2 



ic 7 

6 



ic 8 

2 



(25) 
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Similarly, the decay amplitude for B s — > can be written as: 



M(p u K") = F eK 

+M eK 



+ -C 2 



6 



1 3 15 

— Ca — Ca ~\~ — C*7 -|- — Ca ~\- —Cq ~\- Cin 

3 2 2 3 



/i 



3 13 



+^5 



-c 5 + ~c 7 



fi + M aA - 



/i 



•C3 + 



/a 



+M, 



Pi 



-6 ( -C 5 + -C 7 



/l + Fc 



a A' 



6 --C 3 - C 4 + -C 9 + -C 



10 



<t 



-C5 — C 6 + -C7 + -C 8 
3 o 2 



1; 



fl 

(26) 



where /1 = l/y/2. 

For the decay amplitude of B s — > cjF" decay, one can obtain its decay amplitude from 
Eq. (J2HJ by replacing the vector meson p with u>, i.e, 



(27) 



Note that we have considered the difference in the quark components for the two scalar 
mesons K + and K° in the analytic expressions. We denote the corresponding amplitudes 
for B s -> u/K° decay as F^, M^, M^ 1 , M' aK , M' a %, F' aK and F^ 2 , but do not show 
explicit expressions of these amplitudes here for the sake of simplicity. The total amplitude 
finally can be written as 



MiuoK 



& ( Ci + -c 2 \h-it\ lc 3 + ~c 



Kk 



+2C 5 + jjc 6 + l -C 7 + ic 8 + ^C, - ~C 10 ) h 



+M, 



eK 



-Kk 1 it ( C 5 - -C 7 J h - Kk it ( £C 8 + 



6 (c 5 - \C,\ h - Kk it (Cs - \c^j h 



1 



-c l0 f 



:Cio fi 



(2f 



where /1 = 1 / \/2 . 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



A. Input parameters and wave functions 



In the numerical calculations we use the following input parameters 



A 



MS 

K 




T 



160MeV, 



rn 



250MeV, f p = 205MeV, f\ 
1.6GeV, f Ba = 236MeV, f K = 160MeV, 



m w = 0.782GeV, f u = 200MeV, /J = 160MeV, 

m p = 0.770GeV, M Bs = 5.37GeV, M w = 80.42GeV. (29) 

The central values of the CKM matrix elements to be used in numerical calculations are 



IK 



ud | 



0.9745, \V, 



ub I 



0.0038, \Vt 



tb\ 



\Vtd\ = 0.0083. 



For the B s meson wave function, we adopt the model 
4> Bs {x,b) = N Bs x 2 (l - x) 2 exp 



2ul 2 y 1 



(30) 



(31) 



where ujb is a free parameter and we take Ub = 0.5 ± 0.05 GeV in numerical calculations, 
and N Ba = 65.332 is the normalization factor for uj^ = 0.5. This is the same wave function 
as used in Refs. 0, El Hi . 

For the light meson wave function, we neglect the b dependant part, which is not 
important in numerical analysis. We choose the wave function of p and lu meson as given 
in Ref. ET 



J _p(u) 

2VQ 



f pH x(l-x) l + 0.18C 2 3/2 (2a;-l 



{3(2x - l) 2 + 0.3(2x - l) 2 [5(2x - l) 2 - 3] 



+0.21(3 - 30(2.x - l) 2 + 35(2a; - l) 4 ]} , 
0*0 (*) = - ^) [l + 0.76(10a; 2 -10x + l 

The Gegenbauer polynomial is defined by 



cTit) = $ (5t 2 - 1) 



(32) 

(33) 
(34) 

(35) 



For K meson's wave functions, 0^-, 0^- and 4>k describe the axial vector, pseudoscalar 
and tensor components respectively. We utilize the result from the light-cone sum rule j^] 
including twist-3 contribution: 

(/)£(x) = —f=f K x{l -x) [l + 0.15t + 0.405 (5t 2 - l)] , 



4>k{x) 
f K {x) 



V2N C 
1 

1 



K 



1 + 0.106 (3t 2 - 1) 



0.148 



(3 - 30t 2 + 35t 4 ) 



f K t [1 + 0.1581 (5t 2 -3)] 



(36) 
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whose coefficients correspond to tuq =1.6 GeV and t = 1 — 2x. The parameter tuq is 
defined as: 



m$ = ? ™* , • (37) 

For the pseudoscalar K meson we have assumed that the wave function of K° is the 
same as the wave function of K~, since one can assume that the mass of the u quark is 
equivalent to the mass of the d quark due to the isospin symmetry. 

We include full expression of twist— 3 wave functions for light mesons. The twist— 3 
wave functions are also adopted from QCD sum rule calculations ji^. Using the above 
chosen wave functions and the relevant input parameters, we find the numerical value of 
the corresponding form factor at zero momentum transfer from Eqs. f|14|) : 



F B ^ K (q 2 = 0) = ^ = 0.2761K , (3* 

m B s 

for cjft = 0.50 ± 0.05 GeV, which agrees well with the value as given in Refs. BUS. 



B. Branching ratios 



For B s — > p{u))K decays, the decay amplitudes as given in Eqs. (|23j) . (|2l)j) and 
be rewritten as 



where a = arg 



M = V* b V ud T - V&VttP = V* h V ud T [I + ze % ^] 



can 



(39) 



V u dV* b 



is the weak phase (one of the three CKM angles), and S is the 

relative strong phase between the tree ( "T" ) and penguin ("P" ) amplitude, while the 
term "z" describes the ratio of penguin to tree contributions and is defined as 



(40) 



The ratio z and the strong phase S can be calculated in the pQCD approach. For B s 
p + K~ , p°K and uK decays, we find numerically that 



v: b v td 




p 


V ub V ud 




T 



z(p + K-) 


= 0.10, 


5(p + K-) 


= +138° 




= 1.5, 


5(p K°) = 


= +77°, 




= 1.9, 


5(ujK°) = 


+273°. 



(41) 
(42) 
(43) 



For B s — > pK and B s — > coK decays, the "Tree" diagram contribution is suppressed by a 
factor of Ci + C2/3 ~ 0.1. Thus the penguin diagram contribution is comparable with the 
tree contribution. In our calculation, the only input parameters are wave functions, which 
stand for the non-perturbative contributions. Up to now, no exact solution is made for 
them. So the main uncertainty in pQCD approach comes from these wave functions. In 
this paper, we choose the light cone wave functions which are obtained from QCD Sum 
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Rules |28l . l29j . For heavy B and B s mesons, its wave function is still under discussion 
using different approaches. In this paper, we find the branching ratio of B s — ► p(u>)K is 
sensitive to the wave function parameter ujf,. So the main errors of the ratio z and the 
strong phase S are induced by the uncertainty of Ub = 0.5 ± 0.05 GeV. We just use the 
central values of z and 5 in the following numerical calculations, unless stated explicitly. 

From Eq. (|3*9*|) , it is easy to write the decay amplitude M. for the corresponding charge 
conjugated decay mode. And the CP-averaged branching ratio for considered decays is 
generally defined as 

Br = (\M\ 2 + \M\ 2 )/2 = KbV* d T\ 2 [1 + 22 cos a cos 5 + z 2 ] , (44) 

where the ratio z and the strong phase 5 have been defined in Eqs.([39j) and (|40|). 

Using the wave functions and the input parameters as specified in previous sections, it 
is straightforward to calculate the CP averaged branching ratios for the three considered 
decays. The theoretical predictions in the pQCD approach for the branching ratios of the 



decays under consideration are the following: 

Br(B s -> prtjfF) = [247l$^)3i(«)] x lO" 6 , (45) 

Br(B s - p°K°) = ± 0.1(a)] x 10~ 7 , (46) 

Br(B s ^uoK°) = [l.7lg;|(w 6 ) ± 0.02(a)] x 10~ 7 , (47) 



where the major errors are induced by the uncertainty of uob = 0.5 ± 0.05 GeV, and 
a = 100° ± 20°, respectively. It is easy to see that the above pQCD predictions for the 
branching ratios are consistent with those predicted in the QCD factorization approach 

Br(B s -> p ± K T ) = [24.5i}|;jj] x 10~ 6 , (48) 
Br(B s ^ p°K°) = [6.1™] x 10~ 7 , (49) 
Br(B s ^LuK°) = [5.1+^] x 10" 7 , (50) 

within one standard deviation, although the central values of pQCD predictions for B s — > 
(p°,u)K decays are much smaller than those in QCD factorization approach. 

In Figs. El and El we show the a-dependence of the pQCD predictions for the CP 
averaged branching ratios of B s — ► p^K T , p°K and luK decays for a = [0°, 180°] and 
u b = 0.5 ± 0.05 GeV. 

It is worth stressing that the theoretical predictions in the pQCD approach have large 
theoretical errors induced by the large uncertainties of parameter and CKM angle a. 
From Figs. 121 and El we observe that the pQCD predictions are sensitive to the variations 
of u>b- This sensitive dependence could be fixed by the well measured B s —>■ K form factors 
from the semi-leptonic B s decays as expected in LHCb experiment. Other uncertainties 
in our calculation include the next-to-leading order a s QCD corrections and higher twist 
contributions, which need more complicated calculations. The parameter w 1.6 GeV 
characterizes the relative size of twist 3 contribution to twist 2 contribution. Because of 
the chiral enhancement of raff , the twist 3 contribution become comparable in size with 
the twist 2 contribution. 
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FIG. 2: The branching ratio (in unit of 10~ 6 ) of B s -> p±K T decay for uj h = 0.45 GeV (dotted 
curve), 0.50 GeV (solid curve) and 0.55 GeV (dashed curve), as a function of the CKM angle a. 



C. CP- violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries of B s — > p{uj)K decays 
in pQCD approach. For B s — ► p + K~ decay, the direct CP- violating asymmetry Acp can 
be defined as: 

A dir = \M\ 2 -\M\ 2 = ~ 2z sm a sin 5 
CP \M\ 2 + \M\ 2 1 + 2zcosacos5 + z 2 ' { ' 

where the ratio z and the strong phase 5 have been defined in previous subsection and 
are calculable in pQCD approach. 

Using the central values of z and 6 as given in Eqs. (j4"Tj) , (}4"2j) and (HHJ), it is easy to 
calculate the CP-violating asymmetries: 

Af P (B s - p±K*) = (-12.5115(^)^(0)) x lCr 2 . (52) 

Here two major errors are induced by u>b = 0.50 ± 0.05 GeV and a = 100° ± 20°. The 
pQCD prediction in Eq. (J52)) is also consistent with the prediction in QCDF approach 
|I3 within one standard deviation: Af P {B s -> p ± K T ) = (-1.5 ± 12.2) x 10" 2 . 

In Fig. HJ we show the a— dependence of the direct CP- violating asymmetries A^p 
for B s — > p^K T decay. The possible theoretical errors induced by the uncertainties 
of other input parameters are usually not large, since both z and 5 are stable against 
the variations of them. Uncertainties not included here are the next-to-leading order 
contributions, which may affect the CP asymmetry strongly [3l]|. For B s — > p ± K T decay, 
a large CP asymmetry at 10% level plus large branching ratios at 10~ 5 level are clearly 
measurable in the forthcoming LHCb experiments. 
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FIG. 3: The branching ratio (in unit of 10~ 7 ) of B s — ► p°K and ujK decays for uJb = 0.45 GeV 
(dotted curve), 0.50 GeV(solid curve) and 0.55 GeV(dashed curve), as a function of the CKM 
angle a. 
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FIG. 4: The a-dependence of the direct CP asymmetries of B s — > p decay for uJb = 0.45 
GeV (dotted curve), 0.50 GeV (solid curve) and 0.55 GeV (dashed curve). 



For the pure neutral decays B s — > p°K and uoK , there are both direct and mixing- 
induced CP violation. Using Eq. (joTj) . we find that 

Af P (B s - P°K°) = (-91.9^(u*)3£(a)) x 10~ 2 , 

A d S r P (B s - coif) = (+81.2^; 7 (^):^(a)) x 1(T 2 , (53) 

for u>b = 0.50 ± 0.05 GeV and a = 100° ± 20°. In Fig. |3J we show the a— dependence of 
the direct CP-violating asymmetry for B s — > p°K and B s —>■ uK decays. 
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FIG. 5: The a-dependence of the direct CP asymmetries of B s — ► p°~K° and ujK°) decays for 
= 0.45 GeV (dotted curve), 0.50 GeV (solid curve) and 0.55 GeV (dashed curve). 

The pQCD predictions as given in Eq. are quite different from those obtained by 
using the QCDF approach as given in Refs.jlfjJ : 



A d S r P (B s - P°TC) = (24.71^1) x lCT 2 , 

Afp{B s -> uK°) = (-43.9^;1) x 10" 2 . (54) 

The reason is the great difference in the source of strong phases in two facterization 
approaches. In QCDF approach, the strong phase mainly comes from the perturbative 
charm quark loop diagram, which is a s suppressed But the strong phase in pQCD 
approach comes mainly from non-factorizable and annihilation type diagrams (see figures 
1(c) ~1 (h) ). The sign of the direct CP asymmetry is also different for these two 
approaches for the latter two neutral decays. The forthcoming LHC-b experiments can 
make a test for these two methods. 

Following Ref . 3 > the mixing-induced CP asymmetry for B s — > p !^ and uiK° decays, 
can be defined as 

A mix -2Im(X CP ) sin27 + 2#e(x) sin 7 
CP l + |A C p| 2 l + |x| 2 + 2/?e(x)cos 7 ' { } 

where x = \y Cb Jf* | f^p? an d the angle 7 is one of the three CKM angles. Numerically, the 
pQCD predictions for the mixing induced CP asymmetry are 

A&(B. - p°K°) = (-37^(^)^(7)) x 10- 2 , 
A^(B S - ujK°) = (-40 ± 11(^)^(7)) x 10-2, (56) 

for Lib = 0.50 ± 0.05 GeV and 7 = 60° ± 20°. The 7-dependence of the mixing-induced 
CP asymmetries are shown in Fig. |H1 



16 




30 60 90 120 150 180 



y (degree) 

FIG. 6: The 7-dependence of the mixing-induced CP violating asymmetry of B s — > p°K° (solid 
curve) and B s — ► ujK (dashed curve) for fixed ujb = 0.50 GeV. 



V. SUMMARY 

In this paper, we calculate the branching ratios and CP-violating asymmetries of B s — > 
p ± K If , p°K and cu-K" decays in the pQCD factorization approach. 

Besides the usual factorizable diagrams, the non-factorizable and annihilation diagrams 
as shown in Fig. ^ are also calculated analytically. The non-factorizable and annihilation 
contributions provide the necessary strong phase required by a non-zero CP-violating 
asymmetry for the considered decays. 

From our calculations and phenomenological analysis, we found the following results: 

• From analytical calculations, the form factor F Bs ^ K (0) can be extracted. The 
pQCD prediction is F B ^ K (0) = 0.276^40 for = 0.50 ±0.05 GeV, which agrees 
well with the value as given in Ref. Ill and the result obtained from the QCD sum 



rule calculations [3 

For the CP-averaged branching ratios of the three considered decay modes, the 
pQCD predictions are 

Br(B s -> p ± K T ) « 24.7 x 10" 6 , 
Br(B s -»• p°K°) ss 1.2 x 10" 7 , 

Br(B s -> uSlf) w 1.7 x 10" 7 . (57) 

The theoretical uncertainties are around thirty to fifty percent. The leading pQCD 
predictions for the branching ratios are also consistent with those obtained by em- 
ploying the QCD factorization approach. 

For the CP-violating asymmetries, the theoretical predictions in pQCD approach 
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are 

\dir j jz> , + 



0. 



^(5 S -> p 7r) « -92%, (B a -> p°lT) « -36%, 



-> cjTT) « +81%, -4^?(5 S -> wiT) « -40%, (58) 
for a ~ 100°, 7 ~ 60°. Of course, the theoretical errors here are still large. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function h^s appeared in the expressions of the decay amplitudes in 
Sec. 11111 coming from the Fourier transformations of function H^ ' , 

h e {xi, x 3 , h, b 3 ) = K {^xiX2,m Bs bi) [9(bi - b 3 )K {yfxim Bs bi) h (V%3 m B s h) 

+9(b 3 - b^Ko {^/x~ 3 m Bs b 3 ) I (y/x^msb^] S t {x 3 ), (Al) 

h a (x 2 , x 3 , b 2 , h) = K {i^x 2 x 3 m Ba b 2 ) [9{b 3 - 62)^0 (iV^3 m B s h) h {i\fx~lm B3 b 2 ) 

+9(b 2 - b 3 )K (iy^m B b 2 ) I (iy/x^m Bs b 3 )] S t (x 3 ), (A2) 

h f (xi,x 2 ,x 3 ,b 1 ,b 2 ) = ^9(b 2 - b 1 )l (M Bs ^xix 3 bi)K (M Bs y/x 1 x 3 b 2 ) 

> ( K (M Bs F {1) b 2 ), for ifD>0\ 
+ (h ~b 2 )j. ^ fH (i) (Mss ^Tj for Fh<Q ), (A3) 

h 1 f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = S^9(bi - b 2 )K (^^/x^b 1 M Bs )l (i^x^x 3 'b 2 M Bs ) 

} ( K (M Bs F( 2 )6i), for Ff 2) > 0\ 
+ (h~b 2) y Uh^J^), for F* 2) <oh ( A4 > 



h f {x u x 2l x 3l b Xl b 2 ) = < 9{b x - b 2 )K Q (i^x 2 x 3 b l M Bs )l (i^x 2 x 3 b 2 M j 



B s 



K (M Bs F (3) &i), for F ( 2 3) >0 N 
+ (61 ~ fe) H f^\M Bs M;\b 1 ), for Ffo <oj ' (A5) 
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where Jo is the Bessel function, K and Jo are modified Bessel functions K (—ix) = 

-(n/2)Y (x) + i(ir/2)J (x); H§\z) is the Hankel function ,E. { q\z) = J (z) + i Y (z), and 
F(j)'s are defined by 



1 } 

F(2) 
F (3) 



(xi - x 2 )x 3 , 

xi + x 2 + a; 3 - x x x 3 - x 2 x 3 



The threshold resummation form factor S t (xi) is adopted from Ref.j2 

2 1+2 T(3/2 + c 



VxT(l + c) 



(A6) 



(AT) 



where the parameter c = 0.3. This function is normalized to unity. More information 
about the threshold resummation can be found in reference [l9| . 
The Sudakov factors used in the text are defined as 



s ab (t) 



S cd (t) 



(xim B jV2,bij + s {x 3 m B jV2,hj + s ((1 - x 3 )m B jV2, hj 
1 

Ji 



In Ht / A l I In Ht/A) 



-ln^A) -ln(6 3 A) 
s (xim Bs /V2,bij + s [X2m B jV2,b 2 j + s (jl - x 2 )m Bs /V2,b 

+s (x 3 m Bs /V2, b^j +s((l- x 3 )m Ba /V2, b^j 



1 



2 In 



ln(t/A) 



In 



ln(t/A) 



ln(6iA) ' "~-ln(fo 2 A) 
S ef (t) = s (xtmsj^, bij +s (x 2 m B jV2, b 2 J +s ((1 - x 2 )m Bs /V2~, k 

+s (x 3 m Bs /V2, 6 2 ) + s f (1 - x 3 )m Bs /V2, b : 
1 



Pi 



lnJm + 21n. ln ^ 



-ln(6iA) 



- ln(6 2 A) 



S gh (t) 



= s (x 2 rn B J\/2,b 2 j + s \x 3 m Bs - x 2 )m Bs /V2, 6 ; 

+ sf(l- x 3 )m Ba /V2,k 



1 



]n ln(t/A) + ^ ln(t/A) 



ln(6 3 A) 



ln(M0 



where the function S Bs , S p (^\, Sk used in the amplitudes are defined as: 
S Bs (t) = s(x 1 P+,b 1 ) + 2 f ^(a s (fl)), 

* dp 



S pH (t) = s(x 2 P 2 + , b 2 ) + s ((1 - x 2 )P+, 62) + 2 / ^7 (a fl (/x)) 

•/i/6 2 ^ 

S^(t) = S (x 3 P3-,&3) + s((l-X3)P 3 -,63)+2 / ^7 («.(£)) ■ 



(AS 



(A9) 



(A10) 



(All) 



(A12) 
(A13) 
(A14) 
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where the so called Sudakov factor s(Q, b) resulting from the resummation of double 
logarithms is given as (32! ]: 



*{Q,b)= I ^ 



In 



A(a(fi)) + B(a s (il)) 



(A15) 



with 



A -- 
B 





"67 7T 2 


c F ^ + 


9 3 


71 




e 2 1E -l\ 






3 7T V 





IE 



10 2 _ . / e 
— n; H — on in — 
27 ' 3 V 2 



a. 



71 



(A16) 



Here 7^ = 0.57722 • • • is the Euler constant, n/ is the active quark flavor number. 
The hard scale t^s in the above equations are chosen as: 



max( v /x^m Bs , 1/63) , 
max( v /x7m Bs , l/6 l3 1/6 3 ) , 
max( v /i^m J ; s , l/6 2 , 1/&3) , 
max( v /a^m Bs , I/&2, I/&3) > 
£/ = max( v /xiX 3 m jBs , ^x 2 x 3 m Bs , 1/bi, l/b 2 ) , 
tj. = max(^x 2 x 3 m Bs , I/61, 1/62) , 

= max( v / xi + x 2 + £ 3 - - x 2 x 3 m Bs , y/x 2 x 3 m Bs , l/h, l/b 2 ) . (A17) 

They are given as the maximum energy scale appearing in each diagram to kill the large 
logarithmic radiative corrections. 
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